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Abstract 

We investigate to what extent the temperature dependence of the nuclear symme- 
try energy can affect the neutronization of the stehar core prior to neutrino trapping 
during gravitational collapse. To this end, we implement a one-zone simulation to 
follow the collapse until /3-equilibrium is reached and the lepton fraction remains 
constant. Since the strength of electron capture on the neutron-rich nuclei associated 
to the supernova scenario is still an open issue, we keep it as a free parameter. We 
find that the temperature dependence of the symmetry energy consistently yields 
a small reduction of deleptonization, which corresponds to a systematic effect on 
the shock wave energetics: the gain in dissociation energy of the shock has a small 
yet non-negligible value of about 0.4 foe (1 foe = 10^^ erg) and this result is al- 
most independent from the strength of nuclear electron capture. The presence of 
such a systematic effect and its robustness under changes of the parameters of the 
one-zone model are significant enough to justify further investigations with detailed 
numerical simulations of supernova explosions. 
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1 Introduction 



Weak interaction processes are naturally associated to core-collapse super- 
novae and more generally to compact stars. Indeed, on the one hand the 
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increasing density in the collapsing core of a massive star continuously shifts 
the /^-equilibrium conditions and thus drives electron capture (first mostly 
on exotic nuclei and then on free, unbound protons) all the way to an almost 
completely deleptonized equilibrium state, the final neutron star. On the other 
hand, the tremendous densities and temperatures obtained through the grav- 
itational compression allow the neutrinos, produced both thermally and in 
such weak processes, to interact significantly with matter. The neutrinos dif- 
fuse, rather than stream freely, out of the collapsing core and such neutrino 
transport produces unique physical scenarios, like neutrino trapping half-way 
in the collapse or the shock-wave revival powered by neutrino cooling of the 
proto-neutron star. We refer to the review by Hans Bethe [1] for a masterful 
physical discussion of these phenomena and of their relevance to supernova 
explosions. 

In a previous paper [2] (from now on Paper I; see also Ref. [3] for further de- 
tails), we discussed a particular issue related to electron capture in collapsing 
stellar cores. First, we studied the temperature dependence of the nucleon ef- 
fective mass, m*, in the nuclei ^^Mo, ^"^Zn and ^^Ni as due to the coupling of the 
mean field single-particle levels to the collective surface vibrations of the nu- 
cleus, calculated in the quasi-particle random phase approximation (QRPA). 
Then, we observed that the decrease with temperature obtained for m* in the 
range < T < 2 MeV induces a corresponding increase of the nuclear sym- 
metry energy (per nucleon), -Egym, and in analogy to the results of the Fermi 
gas model we argued that the temperature dependence of -Esym can be fitted 
by a simple analytical expression. Finally, we investigated the implications of 
such a temperature dependence on the gravitational collapse of the core of 
massive stars. We did this in a one-zone (uniform mean density) model, an 
approach which incorporates the important physics but is easy to implement 
and which, in the past, has proven effective to make a preliminary study of the 
core deleptonization during the infall epoch before core bounce, when the col- 
lapse is still homologous [H[5ll6] . In Paper I, electron capture was implemented 
as in the classic work by Bethe et al. (BBAL) [7j, but with the strength of cap- 
ture on nuclei quenched by a factor 7^ = 0.1 to account for the Pauli blocking 
of Gamow- Teller (GT) transitions in neutron-rich heavy nuclei discussed in 
Refs. [5115] . The collapse simulation showed that the temperature dependence 
of the symmetry energy yields a lower rate of neutronization along the col- 
lapse, as expected for larger values of -Egym, and thence a higher value for the 
electron fraction, Y^, at neutrino trapping density. This can be conveniently 
quantified in terms of the associated gain in dissociation ener gy of the shock, 
^T-^diss5 a quantity which gives a more direct physical insighl[^. The results 



A larger lepton fraction after trapping corresponds to a larger homologous core 
so that the shock wave, which forms at its edge after core bounce, will have less 
material to traverse before getting out of the iron core and thence it will dissipate 
less energy in the photo-dissociation of tightly bound nuclei [1]. 
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of Paper I correspond to an energy gain S^Ediss ~ 0.5 — 0.6 foe (1 foe = 10^^ 
erg), a non- negligible amount when one considers that the explosion energy 
(kinetic energy of the ejecta) of SN 1987A was observed to be ~ 1 foe [S]. 

After Paper I was published, two more investigations of the temperature de- 
pendence of the nuclear symmetry energy, due to the presence of neutron-rich 
heavy nuclei in the collapsing core, have been presented p^lITT] : both are 
based on shell model Monte Carlo (SMMC) calculations, the model of choice 
to take into account nuclear correlations beyond those treated at the QRPA 
level. More recently, the virial EOS of hot (T > 2 MeV) low-density nuclear 
matter composed of protons, neutrons and alpha particles was obtained [12]; 
such a scenario is representative of matter around the neutrinosphere after 
core bounce and shock formation, during the contraction phase of the proto- 
neutron star and the associated loss of gravitational energy through neutrino 
emission. The results found in Ref. [12] for the temperature dependence of 
the symmetry energy follow from the presence of these alpha clusters. Such 
findings, however, cannot be applied to the pre-bounce scenario studied here; 
indeed, the temperature of the infalling matter reaches 2 MeV only long af- 
ter /3-equilibrium is achieved and the composition is always dominated by 
neutron-rich heavy clusters in a sea of unbound neutrons, while the fraction 
of alphas is negligible at these low temperatures. The T-dependence of -Egym, 



derived in Refs. [2|10|llj and under discussion here, follows from the presence 

of these exotic nuclear species. 

In Ref. [To], several isobaric pairs with mass numbers in the range A = 54 — 64 
were studied. Although the results obtained for the nuclei ^^Zn and ^^Ni were 
in agreement with those of Paper I, in their conclusions the authors claimed 
to "find no systematic temperature dependence of the symmetry energy coef- 
ficient!!], ^sym, for T' < 1 MeV. This contradicts a recent suggestion that fegym 
increases by 2.5 MeV at this temperature" [10]. An improved SMMC calcula- 
tion, however, was presented several year later in Ref. [IT], where some known 
problems of the previous paper (small model space and g-extrapolation pro- 
cedure to circumvent the notorious sign problem of SMMC) had been fixed. 
Nine isobaric pairs with A = 56 — 66 were analyzed and this time the authors 
concluded that their "SMMC studies are consistent with an increase of the 
symmetry energy with temperature, supporting the argumentation of Donati 
et al." [11]. Indeed, upon averaging over the various pairs, they found a vari- 
ation Sbgym = (6.2 ± 1.8)% in the temperature interval T = 0.33 — 1.23 MeV, 
which is in reasonable agreement with the QRPA results of Paper I, namely 
an increase of the symmetry energy of ~ 8% in the interval T = — 1 MeV. 



^ The coefficient 6sym(T') is the same as the coefficient s{T) of Paper I, and it is 
related to Ssym(^) by the standard expression Egym = ^sym(l — 2a;)^, with x = Z/A 
for a nucleus of mass A and charge Z. 
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In the concluding section of Ref. [TT], the authors also quickly discussed pos- 
sible consequences for core-collapse supernovae. They studied the decrease of 
electron capture on nuclei due to the proposed temperature dependence of the 
symmetry energy, by considering the increase of reaction Q-values induced by 
it. For the neutron-rich nuclei expected during collapse (mass number A > 65), 
they adopted new capture rates obtained in the so-called "hybrid" model [13], 
an approximate approach which mixes SMMC and RPA techniques to go 
beyond the independent particle model (IPM) in the calculation of both al- 
lowed and forbidden transitions. These new results show unblocking of the GT 
strength [13], due to configuration mixing by the residual interaction and to 
thermal excitations, which significantly modify the naive single-particle occu- 
pations of the IPM and thus yield capture rates one order of magnitude larger 
than those predicted by the IPM [5]. Proceeding in this way, the authors ob- 
tained changes of electron capture rates due to the T-dependence of fesym that 
"appear to be rather mild so that one does not expect significant changes for 
the collapse trajectory" [11]. Although we agree that no dramatic effect on 
the dynamics of the collapse is to be expected, one should be more cautious in 
dismissing any significant consequence of the T-dependence of -Egym without a 
collapse simulation. Indeed, not only the reaction Q-values (as considered in 
Ref. [21), but also the equation of state of bulk dense matter, the free nucleon 
abundances, the degree of dissociation into a-particles and the nuclear inter- 
nal excitations are affected by changes in the symmetry energy. Moreover, the 
dynamics of the collapse depends in a very non-linear way on the strength of 
nuclear electron captur^, so that mild changes in the rates may still result 
in non-negligible alterations of the overall energetics. 

The purpose of the present article is to investigate with a collapse simulation 
the extent to which the temperature dependence of the nuclear symmetry en- 
ergy, found in Paper I and confirmed in Ref. [11], can affect the deleptonization 
of the collapsing stellar core. We must, of course, take into account the remark- 
able progress made in the SMMC calculations of electron capture rates since 
publication of Paper I. On the one hand, the new values obtained with im- 
proved SMMC techniques for capture on nuclei present in lower-density matter 
{A < 65) [15] have been implemented in modern evolutionary stellar calcula- 
tions yielding new presupernova models [16], which are significantly different 
than those used so far as initial conditions in collapse simulations. On the 
other hand, the unblocked GT strengths found with the hybrid model for the 
neutron-rich nuclei typical of higher- density matter [A > 65) [T3] have been 
used in numerical (1-dimension) collapse simulations, both newtonian and 
relativistic. When compared to the results from the commonly used Bruenn 
parametrization of nuclear electron capture [T7], which quenches capture on 



^ The parameter study of Ref. [14J, for example, shows that each increase of the 
rate of capture by a factor 10 corresponds roughly to the same decrease (~ O.IMq) 
of the mass of the homologous core. 
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heavy nuclei as required by the IPM and thus allows capture on free (unbound) 
protons to dominate some crucial phases of the collapse, the simulations with 
the new rates show significant differences in the dynamics of the shock wave 
and in the neutrino luminosity 



Altogether, the results obtained in the collapse simulations of Paper I have to 
be revisited in four main aspects, all related to electron capture on nuclei: 
i) hj using the approach of BBAL [7] in Paper I, we certainly overestimated 
the effect of the temperature-dependent symmetry energy on the deleptoniza- 
tion. Indeed, the BBAL rates for capture on nuclei are calculated applying 
the Fermi approximation to a shell model description of the GT transition. 
This statistical limit (which actually does not apply to the collapse scenario, 
where the shell structure is still dominant and the nuclear density of states 
is far from thermal involves an integration over the initial proton states 
and this multiplies the final capture rates by a factor containing the nucleon 
effective maso. This linear dependence of the nuclear rates on m* obviously 
amplifies the thermal effects, but it is absent if a more realistic, non-statistical 
description of capture is adopted. 

a) the BBAL rates for electron capture (on both nuclei and free protons) used 
in Paper I were calculated at T = 0, but since we are looking for a small ther- 
mal effect we cannot neglect the influence of the Fermi distribution functions, 
which describe the occupation numbers of initial and final particle states at 
finite temperature |^|f5] . 

Hi) the multiplying factor 7^ = 0.1, introduced in Paper I to account for the 
Pauli blocking of GT transitions, is not anymore realistic according to the new 
results from the hybrid model |13j. These new findings, however, are not yet 
obtained in a consistent SMMC calculation so that, in our opinion, the actual 
strength of nuclear electron capture is still an open issue and the correct value 
of 7^ is not yet pinned down. 

iv) the initial conditions adopted in Paper I for the collapse have to be revis- 
ited, to account for the new results obtained for the presupernova core when 
implementing the improved SMMC capture rates in evolutionary stellar codes 

m- 



In the next section, we describe our model for the gravitational collapse of the 
stellar core and discuss how it takes into proper account all these issues. 



The integration requires the nuclear density of states, which in the Fermi gas 
model is proportional to the nucleon mass [7] . 
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2 Physical model for the collapse 



In order to study the neutronization of matter induced by gravitational col- 
lapse, we develop a one-zone model (sphere of uniform density) along the 
classic approach of Refs. [HISII^ . The model is an improvement over the one 
used in Paper I in two respects: first, the treatment of electron capture is re- 
visited in order to answer the issues i) and ii) previously mentioned; then, the 
trapping of neutrinos is treated more realistically and provides the equilibrium 
lepton fraction after trapping, when the collapse is adiabatic. Moreover, the 
capture strength on nuclei is kept as a free parameter, 7^, as discussed in issue 
in), and the presupernova initial conditions are the improved ones mentioned 
in issue iv). 

We now describe the main features of our collapse model: 

(1) The dynamical evolution of density with time due to gravity decouples 
from the thermodynamical equations for the changes in entropy and lepton 
fractions. Therefore, we can follow the relevant thermodynamical variables 
(entropy, temperature, electron and neutrino fractions, particle abundances, 
nuclear composition) as a function of density, along the so called collapse tra- 
jectories. 

(2) We adopt the equation of state (EOS) for hot dense matter derived in 
BBAL iTj . The ensemble of nuclear species is approximated by a mean heavy 
nucleuiil in a sea of dripped-out free neutrons and (fewer) protons. The frac- 
tions of free nucleons are determined from nuclear statistical equilibrium. The 
symmetry energy appears in the bulk nuclear energy and, as a consequence, in 
the neutron chemical potential, and in the neutron-proton energy differ- 
ence, fi = fin — f^p- These are crucial quantities in determining the free particle 
abundances, the nuclear capture Q-values and the entropy changes due the 
departure from /^-equilibrium of the collapsing core before neutrino trapping. 

(3) Thermal dissociation of nuclei into a-particles and nucleons is also taken 
into account through the Saha equation, but found to have a negligible effect 
on the collapse trajectories. 

(4) Entropy terms are included for the translational degrees of freedom of all 
the particles (mean heavy nucleus, free classical nucleons, relativistic degener- 
ate leptons) as well as for the internal nuclear excitations, treated in the Fermi 
gas approximation . The nuclear excitation energy is proportional to the nu- 
cleon effective mass (see Ref. [T|), which is the quantity whose temperature 
dependence we originally calculated and fitted by an analytical expression in 
Paper I. We find that the corresponding entropy term has a non-negligible 
effect on the collapse trajectories. 



° An ensemble of nuclei is actually present, in nuclear statistical equilibrium under 
strong and electromagnetic interactions. The mean nucleus is the one that minimizes 
the nuclear energy and thus it represents the most abundant nuclear species [5]. 
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(5) Neutrino trapping is set to start at a given trapping density, ptr- The 
typical "standard" value is ptr,io = 43 (pio being the density in units of 10^° 
g cm~^), but we keep it as a model parameter. As long as p < ptr, neutri- 
nos are allowed to stream freely out of the core and the neutrino fraction is 
Fi, = 0. When p > ptr, neutrino diffusion is treated along the lines of Ref. [6|: 
a degenerate sea of neutrinos with 7^ is allowed to build up by the inclu- 
sion of a diffusion term which decreases with density. Moreover, the inverse 
reactions induced by the sea of neutrinos are included in the electron capture 
rates [20], so that weak interactions can reach equilibrium. In this way, com- 
plete neutrino trapping is reached gradually at a density somewhat larger than 
Ptr ; both the total lepton fraction, Yi = + Y^,, and the entropy tend nat- 
urally to constant values, after which the collapse proceeds adiabatically and 
in /3-equilibrium. This is a major improvement over Paper I, where neutrinos 
were always streaming out freely (Y^ = 0), so that equilibrium could never be 
reached and the equilibrium lepton fraction was just the value of the electron 
fraction taken at p = ptr along the collapse trajectory, namely Ye = Ye^ptr)- 

(6) Electron capture is implemented on both free protons and heavy nuclei 
with standard two-level transitions, as fully developed in Ref. [5]; the phase 
space integral is calculated numerically, although its approximation by Fermi 
integral (as in Eq. (1) of Ref. [TS]) turns out to be accurate enough. For this 
kind of transitions, the nuclear capture rate An is a function of density, tem- 
perature and two other quantities: the excitation energy of the nuclear GT 
resonance. An, and the reaction Q-value. The first is taken as a model param- 
eter, while Q = An (we have actually used a regularized expression for the 
GT excitation energy [20] )• We have also multiplied the nuclear strength An by 
a free parameter, 7^. As shown in Ref. [18], the Q-dependence of the capture 
rates obtained with the hybrid model can be reasonably fitted by the two-level 
expression, with An = 2.5 MeV and an appropriate GT matrix element; we 
normalize An so that our expression coincides with Eq. (1) of Ref. [I8j when 




(7) The temperature dependence of the symmetry energy is treated as in 
Paper I pi's], where it was expressed in terms of the T-dependence of the nu- 
cleon effective mass, m* = m*{T), calculated for different nuclei. The results 
for each nucleus were fitted with a formula containing two parameters: the 
value at T = of the so-called cu-mass, m^{0), and the temperature scale of 
this dependence, Tq. The standard average values are m^{0) = 1.7 and Tq = 2 
MeV, but we keep them as model parameters allowed to vary in a meaningful 
physical range (1.4 < m^{0) < 1.8 and 1.9 < Tq < 2.1 MeV [2]), to account 
for their dependence on the nucleus studied. 

Details and equations can be found in Refs. [21[3|lf51l5f7] . In particular, we 
adopted the following differential equations for the collapse trajectories: 
i) the electron fraction evolution, dYe/dp, is Eq. (85) of Fuller [5]. From trap- 
ping density onwards, the neutrino-induced inverse reactions are included as 
in Eqs. (15) and (16) of Ray et al. ^. 
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a) the entropy and temperature evolution, els'/dp and dT/dp, are respectively: 
Eqs. (92) and (93) of Fuller before trappin£3; Eqs. (33) and (35) of Ray et al. 
after trapping. 

in) the neutrino fraction evolution, dFj,/dp, is Eq. (29) of Ray et al. after 
trapping; before trapping we take Yi, = 

The collapse trajectories are determined starting from a set of initial conditions 
on the density, pi, the temperature, Tj, and the electron fraction, Fe.i (until 
trapping density is reached, Yy = 0). According to the improved results of 
Ref. [IH] for the central properties of the presupernova core which evolves 
from a ISM© star (about the size of the progenitor of SN 1987A), we will take 
the initial values pio,i = 0.936, T; = 0.625 MeV and Ye,i = 0.432, which differ 
significantly from those adopted in Paper I. The differential equations are 
then integrated and the collapse trajectories of the different thermodynamical 
quantities are found. In particular, the total lepton fraction Yj = Yi{p) tends to 
a constant value, F/^tr, as the density increases above ptr and neutrino trapping 
is completed. 

In the next section, we discuss our results for the neutronization of the core 
in terms of tr and of quantities related to it. 



3 Results of the collapse simulation 

We first fix the model parameters to their "standard" values (ptr,io = 43, 
^u;(0) = 1-7, Tq = 2 MeV, An = 2.5 MeV) and make a parameter study 
of the core neutronization as a function of the nuclear strength in the range 
< 7^ < 5. We point out that 7^ = corresponds to electron capture 
on free protons only, while 7^ = 5 is very large and probably unrealistic. 
The older, blocked GT rates of Fuller [5] correspond to 7^ = 0.1, while the 
new unblocked rates of Ref. [IB] are associated to 7^ = 1. Improved future 
calculations could change the presently accepted value of the nuclear strength, 
but (barring discovery of past errors or unexpected breakthroughs) we think 
that 0.5 < 7^ < 2 should represent a reasonable physical range. 

For each choice of parameters, we have run the collapse simulation twice: once 
implementing the temperature dependence -Egym = EsymiT) and obtaining 
yj^trly, once setting -Egym = -E'sym(O) and obtaining Vz,trlo- We indicate by 6^ 
the "thermal" variation of a quantity due to the temperature dependence of 
the symmetry energy; for example, the thermal change in equilibrium lepton 
fraction is 5j,YJ,tr = y«,trlT - ^,trlo- 



^ In Eq. (92) of Fuller, the terms AXp in the denominators should each be multiplied 
by a factor Afp. 
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From the collapse trajectories of the lepton fraction, Yi = Yi{p), we can infer 
the density range where the temperature dependence of the symmetry energy 
is most relevant. Regarding the thermal variation of the capture rates, we ob- 
tain similar results as those shown in Figure 7 of Ref . [TI] : the thermal effect 
is maximum at the onset of collapse (pio ~ 1) and decreases with increas- 
ing density; when trapping sets in (pio ~ 40), the difference in the capture 
rates due to EsymiT) is already negligible. However, since this is a cumulative 
and non-linear effect on the lepton fraction, a difference between the collapse 
trajectories, S^Yi{p) = Yi{p)\^ — Yj(p)|Q, gradually builds up as the density 
increases. This is shown in Figured], where the thermal variation of the lepton 
fraction collapse trajectories, 6j.Yi, is given as a function of the density of the 
collapsing core. We notice that, even though the thermal variation actually 
reaches its final equilibrium value, 5^ Yj^tr, only when neutrino trapping is fully 
achieved (around pio ~ 250), the whole effect is seen to build up before trap- 
ping sets in. In particular, it is mostly in the low-density range 1 < pio < 20 
that the temperature dependence of the symmetry energy affects significantly 
the neutronization process. 
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Fig. 1. Thermal variation of the lepton fraction collapse trajectories, S^^Yi, as a 
function of the density of the collapsing core, pio (in units of 10^'' g cm~^). The 
dotted line indicates the trapping density, ptr- The curve corresponds to standard 
parameters of the model (ptr,io = 43, m^{0) = 1.7, Tq = 2 MeV, An = 2.5 MeV). 

Surprisingly, although Figure [1] corresponds to standard parameters of the 
model, the curve 5j,Yi(p) is found to change very little under physically rea- 
sonable variations of the model parameters, as will be discussed later on. 

The results for the equilibrium lepton fractions and the associated thermal 
changes are shown in Table [1] for different values of 7^ . We point out how the 
general magnitude of the equilibrium lepton fraction is a very slowly decreas- 
ing function of 7^. Increasing the strength by a factor ten, from the blocked 
to the unblocked capture rates, decreases the equilibrium lepton fraction by 
~ 13%, which is in reasonable agreement with the ~ 10% change obtained 
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in newtonian one-dimensional simulations [H]. We also notice that the ther- 
mal effect under study systematically reduces the equilibrium neutronization, 
namely Yj^tr is increased by an almost constant value, 5j,yi,tr — 0.006, irrespec- 
tive of the value of the strength 7^. Although small, this effect is not negligible, 
as we will argue in the remaining of this article. 
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Table 1 

Results of the collapse simulation for different values of the strength of nuclear elec- 
tron capture, 7^. We show the equilibrium lepton fractions after trapping obtained 
without, ll,trlo) and with, l/^trlri the temperature dependence of the symmetry en- 
ergy; the thermal change of the equilibrium lepton fraction, 5j,Yi^i^; the correspond- 
ing gain in dissociation energy of the shock, (^^ii^diss (in foe). The model parameters 
are the standard ones (/9tr,io = 43, m^(0) = 1.7, Tq = 2 MeV, An = 2.5 MeV). 

In order to determine the relevance of our results to supernova explosions, we 
need a quantity with a more direct physical meaning and which can be com- 
pared to relevant observables. As in Paper I, we use the gain in shock dissoci- 
ation energy which is defined as 5j,i?diss = 98 [(yz,tr|T)^ ~ (^«,tr|o)^] = 98 6j.Yi\^ 
(in foe). Although based on a schematic model for the shock formation and 
propagation [21j, this expression provides a reasonable order of magnitude 
estimate of S^E^iss- In a similar fashion, one could consider the change in ini- 
tial (i.e. post-bounce) shock energy, ^^ii^shock, which also foUows from changes 
in the equilibrium lepton fractions affecting the size of the homologous core. 
In the schematic approach of Ref. [Sj, however, the expression for the initial 
shock energy, -Eshock = -Eshock(yi,tr), has a maximum for YJ^tr = jf^i = 0.3323. 
Since the equilibrium lepton fractions corresponding to 7^ = 1 are close to 
this extremum (cf. Table [1]), the thermal effect (5j,-Eshock turns out to be quite 
small (~ 10~^ foe); we will not consider it in the following. 

Since ^^YJ^tr is small, the thermal gain in dissociation energy can be written 
as 5j,i?diss — 196 Y;,trlo ^ ^T^i,tr- This shows that in general 5j,i?diss depends 
on S^Yi^tr, but its magnitude is fixed by the final neutronization reached by 
matter, yj,trlo5 which is determined by the nuclear capture strength 7^. In the 
last column of Table [H we show the results for the gain in dissociation energy. 
For standard parameters and 7^ = 1, we find (5y£'diss = 0.39 foe. Moreover, 
since S^Yi^tr is constant, the gain in dissociation energy has the same very slow 
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dependence on the strength parameter as the equihbrium lepton fraction. This 
is well seen in Figure [21 where ^^E'diss is given as a function of 7^. The points 
are the results of the collapse simulation, while the line in the log-log graph 
represent a power-law best fit, with a very small exponent m = —0.065. In 
the physical meaningful range for the strength (0.5 ^ 7^ ^ 2), the gain in 
dissociation energy varies only by ±4%, in the interval 6^Ediss ~ 0.38 — 0.41 
foe. 




0.1 1 ■ • — ■ — ■ • — ■ — 

0.1 0.5 1 5 10 

2 

y 

Fig. 2. Gain in dissociation energy of the shock, ^yE'diss (in foe), as a function of 
the strength of nuclear electron capture, 7^. The calculated points correspond to 
standard parameters of the model (ptr,io = 43, muj{0) = 1.7, Tq = 2 MeV, An = 2.5 
MeV). The line represents a power-law fit, with exponent m = —0.065. 

Although the previous discussion indicate a quite stable value ^^i^diss ~ 0.4 
foe, we want to study the robustness of such a result under reasonable varia- 
tions of the model parameters, compatible with present theoretical uncertain- 
ties about the values of ptr,io, An, muj{0) and Tq. In Table O we show S^E^iss 
for different values of 7^: each column represents the case in which only one of 
the parameters, ptr,io or An, is changed from its standard value to the value 
indicated. In Figure [3l instead, we show a contour plot for Sr^E^iss (in foe) as 
a function of the two parameters mj^(O) and Tq, the other ones being fixed at 
their standard values. The solid level lines are for 7^ = 1, the dotted ones for 
7^ = 0.1 and the shaded area indicates the physically meaningful range found 
in Paper I for the thermal parameters of the symmetry energy. 

The results of Table |2] and Figure [3] show that, under reasonable variations 
of the model parameters, the gain in dissociation energy of the shock changes 
only by about ±10%, in the range Sr^E^iss ~ 0.35 — 0.45 foe. This proves the 
robustness of our conclusions: the temperature dependence of the symmetry 
energy yields a systematic energy gain (less dissipation of shock energy), whose 
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Table 2 

Dependence of the results from the parameters of the model. We show the gain in 
dissociation energy of the shock, (5j,-Ediss (in foe), for different values of the strength 
of nuclear electron capture, 7^. In each column we change only the value of one 
parameter, either the excitation energy of the GT resonance, An (in MeV), or the 
trapping density, /Otr.iO! while the other parameters are the standard ones (ptr.io = 
43, m<^(0) = 1.7, To = 2 MeV, An = 2.5 MeV). 
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Fig. 3. Gain in dissociation energy of the shock, 5j,i?diss (level lines labelled in foe), 
as a function of the parameters 771^^(0) and Tq. The solid contour lines correspond to 
7^ = 1, the dotted contour lines to 7^ = 0.1. The other parameters are the standard 
ones (/Otr,io = 43, An = 2.5 MeV). The shaded area shows the physical range found 
in Ref. [2j for the parameters ?n.(^(0) and Tq. 

order of magnitude is ^^-Ediss ~ 0.4 foS In the concluding section, we will 

As in Paper I, we have set the volume and symmetry energy coefficients in the 
BBAL EOS to the values = -16.5 MeV and s(0) = 29.3 MeV respectively. 
Different values give a different energy gain, but do not alter our general conclusions 
(for example, with wq = —16 MeV and s(0) = 31.3 MeV we find S^E^iss ~ 0.3 foe). 
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discuss the relevance of such a result for the physics of supernova explosions. 



4 Conclusion 

In this article we have studied the effect of the temperature dependence of 
the symmetry energy on the neutronization processes occurring during grav- 
itational collapse in a supernova explosion. We have assumed for -Egym the 
T-dependence found in Paper I and later confirmed in Ref. [11], first fixing 
the parameters to their average values m^(0) = 1.7, Tq = 2 MeV, but later 
allowing them to vary in a reasonable physical interval. We have followed the 
collapse with a one-zone model, finding the collapse trajectories of the differ- 
ent thermodynamical variables (temperature, entropy, particle abundances, 
lepton fractions, mean heavy nucleus) and determining the equilibrium lep- 
ton fraction after neutrino trapping, when the collapse becomes adiabatic. We 
have implemented electron capture on both free protons and nuclei with stan- 
dard two-level transitions at finite temperature. However, to account for the 
present theoretical uncertainties concerning electron capture rates in exotic 
nuclei, we have multiplied the nuclear strength recently obtained in Ref. [18] 
by a strength parameter 7^: variations in a range 0.5 ^ 7^ ^ 2 around the 
presently accepted value of 7^ = 1 are not to be ruled out in the future. 
Starting from the improved presupernova initial conditions of Ref. [16], we 
have run the collapse simulation with and without the T-dependence of £'sym 
implemented, thus obtaining the "thermal" change in deleptonization, 5^,1/,^- 
Then, we have studied the significance of this thermal effect in terms of a quan- 
tity with more direct physical meaning, the corresponding gain in dissociation 
energy of the shock, 5yi?diss oc 5,^Y^^^. Finally, we have tested the solidity of 
our results by varying the standard parameters of the model (ptr,io = 43, 
^(^(0) = 1.7, To = 2 MeV, An = 2.5 MeV) within reasonable physical ranges, 
compatible with present theoretical uncertainties. 

The main conclusion of our investigation is that the temperature dependence 
of the symmetry energy systematically reduces the neutronization of the core, 
namely it consistently increases the equilibrium lepton fraction by a small 
constant amount, 5,^Yi^tv — 0.006, irrespective of the value of the strength 
parameter 7^. The corresponding gain in shock dissociation energy, instead, 
decreases with increasing nuclear strength, but very slowly: when 7^ = 1 is 
divided or multiplied by two, 5y£^diss varies only by ±4% around its standard 
parameter value 5yi?diss| 2_i — 0-39 foe. Moreover, significant changes in the 
other one- zone model parameters correspond to a quite small range o^-Z^/diss ^ 
0.35 — 0.45 foe. This confirms the robustness of our results and the presence 
of a systematic gain in shock dissociation energy of order S^E^iss ~ 0.4 foe, 
associated to the temperature dependence of the symmetry energy. 
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Such an effect is obviously not a dramatic one, when one considers that the 
total energy sapped from the shock by photo-dissociation of nuclei is larger by 
almost two orders of magnitude. Indeed, even changing the nuclear strength by 
a factor of ten through the unblocking of GT transitions does not qualitatively 
alter the final outcome of the failed explosion, at least in one- dimensional sim- 
ulations [19j. Actually, recent developments of three-dimensional simulations 
of core-collapse supernovae indicate that the roles of neutrinos, fluid instabil- 
ities, rotation and magnetic fields are probably critical to obtain successful 
explosions [22] ■ However, when compared to the typical kinetic energies of 
a supernova explosion, -R'expb which are imparted by the shock wave to the 
ejecta, a gain in shock energy of Sr^Ediss ~ 0.4 foe is not negligible (for SN 
1987A, observation gave -K'expi ~ 1 foe [9]). Moreover, 6,^Ediss is two orders 
of magnitude larger than the total electromagnetic output p. On general 
grounds, since both the explosion energy -f^expi and the much smaller elec- 
tromagnetic output have small values resulting from differences of very large 
quantities (gravitational energy, initial post-bounce shock energy, neutrino 
losses, nuclear photo-dissociation), it follows that the explosion observables 
can be sensitive to subtle microphysical features. In particular, systematic nu- 
clear effects can be of particular importance, as noted also in the conclusions 
of Ref. [m. 

The numerical results of our one-zone collapse simulation are significant for 
their order of magnitude, not their precise values which are limited by the over- 
simplified zero-dimensional approach. In our opinion, their robustness under 
variations of the model parameters justifies further investigation in detailed 
one-dimensional numerical codes. This is not a straightforward task, since not 
only the reaction Q-values, but the whole EOS describing dense hot mat- 
ter is affected by the T-dependence of -Egym- In particular, this is lengthy to 
implement within the Lattimer and Swesty EOS |23j, currently used in re- 
alistic supernova codes, since changing any of its nuclear input parameters 
(in the present case introducing m^(0) and To to parameterize m* = ui^iT) 
and thence the temperature dependence of the symmetry energy) requires re- 
calculating the phase and Maxwell construction boundaries before generating 
the final table of the EOS. We are presently working along these lines j24j , 
with interesting preliminary results from one- dimensional simulations with 
the simpler BBAL EOS, where the temperature dependence of the symmetry 
energy can be implemented analytically. 
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